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ABSTRACT 

The spatial cluster-galaxy correlation amplitude, B gc , is computed for a set 
of 76 z < 0.3 ultraluminous infrared galaxies (ULIRGs) from the 1-Jy sample. 
The B gc parameter is used to quantify the richness of the environment within 0.5 
Mpc of each ULIRG. We find that the environment of local ULIRGs is similar to 
that of the field with the possible exceptions of a few objects with environmental 
densities typical of clusters with Abell richness classes and 1. No obvious trends 
are seen with redshift, optical spectral type, infrared luminosity, or infrared color 
(■/W/k))- We compare these results with those of local AGNs and QSOs at 
various redshifts. The 1-Jy ULIRGs show a broader range of environments than 
local Seyferts, which are exclusively found in the field. The distribution of ULIRG 
5 9C -values overlaps considerably with that of local QSOs, consistent with the 
scenario where some QSOs go through a ultraluminous infrared phase. However, 
a rigorous statistical analysis of the data indicates that these two samples are not 
drawn from the same parent population. The B gc distribution of QSOs shows 
a distinct tail at high I? 9C -values which is not apparent among the ULIRGs. 
This difference is consistent with the fact that some of the QSOs used for this 
comparison have bigger and more luminous hosts than the 1-Jy ULIRGs. 

Subject headings: galaxies: active - galaxies: clusters: general - quasars: general 

1. Introduction 

Ultraluminous Infrared Galaxies (ULIRGs) are defined as galaxies with Ltr = L(8 — 
1000 yum) > 10 12 L Q (see reviews by Sanders & Mirabel 1996; Lonsdale, Farrah, & Smith 
2006). This luminosity limit is roughly equivalent to the minimum bolometric luminosity 
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of QSOs. At luminosities above 1O 12 L , the space density of ULIRGs in the local universe 
is greater than that of optically selected quasars with similar bolometric luminosities by a 
factor of ~ 1.5. Thus local ULIRGs represent the most common type of ultraluminous galaxy. 
Systematic optical and near-infrared imaging surveys have revealed that local ULIRGs are 
almost always undergoing major mergers (e.g., Surace & Sanders 1999; Surace, Sanders, 
& Evans 2001; Scoville et al. 2000; Veilleux et al. 2002, 2006). Most of the gas and star 
formation (and AGN) activity in these systems are concentrated well within the central kpc 
(e.g., Downes & Solomon 1998; Soifer et al. 2000, 2001). Ground-based optical and near- 
infrared spectroscopic studies of these objects have shown that at least 25% - 30% of them 
show genuine signs of AGN activity (e.g., Kim, Veilleux, & Sanders 1998; Veilleux, Kim, 
& Sanders 1997, 1999; Veilleux, Sanders, & Kim 1999). This fraction increases to ~ 50% 
among the objects with log[LiR/L Q ] > 12.3. These results are compatible with those from 
mid-infrared spectroscopic surveys (e.g., Genzel et al. 1998; Lutz et al. 1998; Lutz, Veilleux, 
k Genzel 1999; Rigopoulou et al. 1999; Tran et al 2001). 

ULIRGs are relevant to a wide range of astronomical issues, including the role played 
by galactic mergers in forming some or all elliptical galaxies (Genzel et al. 2001; Veilleux 
et al. 2002), the efficiency of transport of gas into the central regions of such mergers and 
the subsequent triggering of circumnuclear star formation (e.g., Mihos & Hernquist 1996; 
Barnes 2004), the resulting heating and metal enrichment of the IGM by galactic winds (e.g. 
Rupke, Veilleux, & Sanders 2002, 2005ab; Veilleux, Cecil, & Bland-Hawthorn 2005; Martin 
2005), the potential growth and fueling of supermassive black holes and the possible origin 
of quasars (Sanders et al. 1988). The discovery of z = 1 — 4 submm sources with SCUBA 
(e.g., Smail et al. 1997; Hughes et al. 1998) suggests that ULIRGs are also relevant to the 
dominant source of radiant energy in the universe today. Indeed, integration of the light 
from the SCUBA population shows that it may account for most of the sub mm/far- infrared 
background, as a result of the strong cosmological evolution of these sources (e.g., Chapman 
et al. 2005). Thus, while the present-day ULIRGs provide a relatively small contribution to 
the total present background, their cousins at high z are fundamentally important in this 
regard. 

If ULIRGs are the predecessors of QSOs, one would expect ULIRGs and QSOs to 
live in similar environments. Surprisingly little has been published on the environments of 
local ULIRGs, in stark contrast to the abundant literature on the small- and large-scale 
environments of AGNs and QSOs (e.g., Yee, Green, & Stockman 1986; Yee & Green 1987; 
Ellingson et al. 1991; Hill & Lilly 1991; de Robertis et al. 1998; McLure & Dunlop 2001; 
Wold et al. 2000, 2001; Barr et al. 2003; Miller et al. 2003; Kauffmann et al. 2004; Sochting 
et al. 2004; Wake et al. 2004; Croom et al. 2005; Waskett et al. 2005; Serber, Bahcall, & 
Richards 2006) and the growing literature on the environments of z > 1 ULIRGs (e.g., Blain 
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et al. 2004; Farrah et al. 2004, 2006). To our knowledge, Tacconi et al. (2002) is the only 
published study that has attempted to quantify the environments of local ULIRGs. They 
correlated the positions of local ULIRGs with the catalogs of galaxy clusters and groups 
available in NED and found that none of them are located within a galaxy cluster. The lack 
of comprehensive imaging database at the time prevented them from carrying out a more 
quantitative clustering analysis of these objects. 

The present paper remedies the situation by using the large imaging database of Veilleux 
et al. (2002) to quantify the environment of local ((z) ~ 0.15) ULIRGs from the 1-Jy sample. 
We note that the spectroscopy portion of the Sloan Digital Sky Survey (SDSS) provides 
redshift information for only the bright tail of the galaxy luminosity function at z ~ 0.15, 
so a method that relies solely on the photometric measurements of the galaxies in the field 
surrounding the ULIRG must be used for the present analysis. The properties of the 1- 
Jy sample and imaging dataset are reviewed in §2. In §3, the procedure for deriving the 
environmental richness, B gc , is outlined. Results for our sample are presented in §4. The 
findings of environmental studies for quasars and Seyferts are compared with our results in 
§5. Our conclusions are summarized in §6. We use Hq = 50 km s _1 Mpc -1 , Q m = 1, and Q\ 
= throughout this paper. These values were selected to match those of previous studies 
and facilitate comparisons; they have no effect on our conclusions. 



2. Sample 

The IRAS 1-Jy sample of 118 ULIRGs identified by Kim & Sanders (1998) is the starting 
point of our investigation. The 1-Jy ULIRGs were selected to have high galactic latitude (|6| 
> 30°), 60-/zm flux greater than 1 Jy, 60-/im flux greater than their 12- /xm flux (to exclude 
infrared-bright stars), and ratios of 60-/zm flux to 100- /im flux above 10~ 03 (to favor the 
detection of high-luminosity objects). 

All 1-Jy ULIRGs were imaged at optical (R) and near-infrared (K 7 ) wavelengths using 
the U. of Hawaii 2.2-meter telescope. The present study uses only the R-band images since 
they have a larger field of view (FOV) and are deeper than the K'-band images. The R filter 
at 6400 A was a Kron-Cousins filter. Details of the observations and data reduction can be 
found in Kim, Veilleux, & Sanders (2002). The analysis of these data is presented in Veilleux 
et al. (2002). These data are part of comprehensive imaging and spectroscopic surveys which 
also include a large set of optical and near-infrared spectra of the nuclear sources (Veilleux et 
al. 1999ab and references therein), a growing set of spatially-resolved near-infrared spectra 
to study the gas and stellar kinematics of the hosts (Genzel et al. 2001; Tacconi et al. 2002; 
Dasyra et al. 2006a, 2006b), and mid-infrared spectra from the Infrared Space Observatory 
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(ISO) and the Spitzer Space Telescope (SST; e.g., Genzel et al. 1998; Veilleux et al. 2006b, 
in prep.). This effort is called QUEST: Quasar / ULIRG Evolutionary Study. 

Since the set of data presented in Kim et al. (2002) was compiled from observations 
made over the course of 14 years, a variety of CCDs were used and the FOV sizes and spatial 
resolutions are not uniform. For consistency, we limit the set of data in this paper to the 
images of the 76 objects taken under good photometric conditions with the TEK 2048 x 
2048 CCD. Of these images, 32 (42%) were irrecoverably cropped during an earlier stage of 
data reduction and have a significantly reduced FOV size. The effects of the cropping on 
the results of our analysis are discussed in §4. Table 1 lists the objects in our sample along 
with the FOV size and several other properties of the sources. 

3. Analysis 

In this section we explain the methods that we used to quantify the environment richness 
around each ULIRG. First, we describe the algorithms used to find objects in the field and 
identify them as stars or galaxies. Next, we discuss the formalism applied to calculate the 
environment richness parameter, B gc . The techniques used for our analysis have already 
been described in detail in Yee (1991), Ellingson et al. (1991), Yee & Lopez-Cruz (1999), 
and Gladders & Yee (2005); here we highlight the main steps. 

3.1. Object Identification and Classification 

Object identification was accomplished using the Picture Processing Package (PPP) 
developed by one of us (Yee 1991). This program systematically examines each pixel in 
the image and determines whether it has the potential to be part of an object: a star, a 
galaxy, a cosmic ray, or an artifact of the CCD. After running through a series of tests, 
the PPP object finding program identifies and catalogs the location and peak brightness of 
objects in the image. The algorithms used here are modified versions of that used by Kron 
(1980), which depend on searching for local maxima. They have been shown to be robust 
for object identification in sparse to moderately crowded fields (Yee 1991). The 1-Jy sample 
selection criterion \b\ > 30° avoids extremely crowded fields (and reduces the effects of dust 
extinction on the galaxy counts), which could lead to object misclassification and erroneous 
environment richness measurements. We therefore find that this object finding routine is 
perfectly adequate for all ULIRGs in our sample 

To address the problem of bad pixels or cosmic rays, objects were thrown out automat- 
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ically if a given pixel was five times brighter than those immediately surrounding it. This 
did not always work well because bright bad pixels are sometimes surrounded by other bad 
pixels. So, some misidentified objects were also identified by eye and removed by hand. 

The next step was to run an aperture photometry algorithm on the identified objects 
in each image to determine whether these objects are stars or galaxies. For each object, a 
growth curve was calculated using a series of circular apertures centered on the intensity 
centroid of the object. A reference-star growth curve was created for each quadrant of 
the CCD frame by averaging the growth curves of bright, isolated, and unsaturated stars 
within each quadrant. The growth curves of the other objects were then compared with the 
reference-star growth curve using the classification parameter C2 defined by Yee (1991). In 
essence, C 2 computes the average difference per aperture between the growth curves of the 
objects and the growth curve of the reference star after they have been scaled to match at 
the center and effectively compares the ratio between the fluxes in the center and the outer 
part of an object with that of the reference star. This method has been thoroughly tested 
by Yee (1991); readers interested in knowing more about this classification scheme should 
refer to this paper for detail. 



3.2. Environment Richness Parameter 

We use the parameter B gc to quantify the richness of the environment of ULIRGs. B gc is 
the amplitude of the galaxy-galaxy correlation function calculated for each object of interest 
individually. It was first used by Longair & Seldner (1979) to measure the environment 
of radio galaxies using galaxy counts, and subsequently adopted in most studies of the 
environments of quasars and other active galaxies (e.g., Yee & Green 1984; Ellingson et al. 
1991; de Robertis et al. 1998; McLure k Dunlop 2001; Wold et al. 2000, 2001; Barr et al. 
2003; Waskett et al. 2005), and also used as a quantitative measurement of galaxy cluster 
richness (e.g., Andersen & Owen, 1994; Yee & Lopez-Cruz 1999). Yee & Lopez-Cruz (1999) 
have demonstrated the robustness of the B gc parameter when galaxies are counted to different 
radii and to different depth. Furthermore, measurements of the environmental richness based 
on the photometrically-derived _B 9C -values have been shown to be entirely consistent with 
measurements based on spectroscopic data. This was demonstrated by Yee & Ellingson 
(2003), who used the data from the Canadian Network for Observational Cosmology Cluster 
Redshift Survey (CNOC1) to compare -B 9C -values derived from (1) photometric data with 
background subtraction, and (2) from properly weighted spectroscopy data to account for 
incompleteness. We describe briefly the procedure for deriving B gc below. 

In order to determine the richness of the environment around a ULIRG, we need to 
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count the number of galaxies within a spherical volume with radius, r, from the ULIRG 
of interest. However, we necessarily must begin with a two-dimensional image, which is a 
projection of this volume onto the sky plane. The number of galaxies in a solid angle dfl, at 
an angular distance 9 from the object of interest is given by (Seldner & Peebles 1978) 

N(9)dQ — N g [l + uj(9)]dQ, (1) 

where N g is the average surface density of galaxies and w(9), the angular correlation function, 
can be expressed approximately as a power law, 

w(9) = A gc 9 l -\ (2) 

A gc is a measure of the average enhancement of galaxies in angular area, and 7 ~ 1.77 
empirically. Integrating equation (2) within a circle with radius 9 yields 

- Ntot — Nb gc (3 — 7) 7 _ x , , 

Age = 7- 9 , [6) 

Nbgc * 

where N tot and N bgc are the the total numbers of galaxies and background galaxies, respec- 
tively, within an angular radius of 9. 

Next, the two-dimensional parameters must be translated into three dimensions. The 
angular correlation function w(9) is translated into the spatial correlation function, £(r), 
which describes the number of galaxies in volume element dV at distance r from the object 
of interest. It can be shown that £(r) = -B 9C r" 7 , where 7 has the same value as in equation 
(3) and B gc is the spatial correlation amplitude, a measure of the richness of the environment 
around the galaxy. Longair & Seldner (1979) have shown that 

Ag C n bg (m)D^ 
9C ~ I^(m,z) ■ [) 

The constant I 7 is an integration constant which depends on 7 (Groth & Peebles 1977). 
nb g (m) is the expected count per unit angular area of background galaxies brighter than 
apparent magnitude to, \1/(to, z) is the normalized integrated luminosity function of galaxies 
to apparent magnitude to, at redshift z of the cluster, and D is the angular diameter distance 
to the ULIRG at redshift z. For our calculations, we used 7 = 1.77, J 7 = 3.78, and the 
cosmological parameters Ho = 50 km s _1 Mpc _1 , Q m = 1, and Q\ — to match those of 
previous papers. For \1/ and n b g(m), we use the luminosity function and background counts 
derived from the Red-Sequence Cluster Survey (RCS; e.g., Gladders & Yee 2005). 

The uncertainty on B gc is computed using the formula 



AB 9C _ (N net + 1.3 2 N bg )2 

D gc iv net 
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where N net is the net counts of galaxies over the background counts, Nj, g . This is a conser- 
vatively large error estimate as it includes the expected counting statistics in N net and the 
expected dispersion in background counts. The factor 1.3 2 is included to account approxi- 
mately for the additional fluctuation from the clustered (and hence non-Poissonian) nature 
of the background counts (discussed in detail in Yee, Green, & Stockman 1986). We follow 
the prescription of Yee & Lopez-Cruz (1999) and integrate the luminosity function from 
approximately M R = —25 to + 2 (where « —22.3 for our cosmology) to calculate 
the galaxy counts. This corresponds roughly to R = 15 - 20 for the galaxies in our sample 
((z) 0.15). This range of integration was found by Yee & Lopez-Cruz (1999) to reduce 
the sensitivity to small intrinsic variation of M* and variations in the faint-end slopes of 
the cluster luminosity function. The B gc parameter is computed over a radius r = 500 kpc, 
either directly from the data when FOV > 1 Mpc or extrapolated to this radius when FOV 
< 1 Mpc. This radius is selected to match that of previous studies. The B gc parameter is 
not sensitive to this radius (Yee & Lopez-Cruz 1999; see also §4). 



4. Results 

The spatial correlation amplitude parameter, B gc , was computed for each of the 76 
ULIRGs in our sample. The B 9C -values and associated 1-a uncertainties are listed for each 
object in Table 1. The average (median) value of B gc and 1-a scatter around the mean 
for our sample of 76 ULIRGs is (B gc ) = 35 ± 198 Mpc 1 ' 77 (—3 Mpc 1 ' 77 ). For comparison, 
the -B sc -values of field galaxies and clusters of Abell richness class (ARC) 0-4 are ~ 67.5, 
600 ± 200, 1000 ± 200, 1400 ± 200, 1800 ± 200, and 2200 ± 200 Mpc 1 - 77 , respectively 
(the field _B 9C -value is from Davis & Peebles 1983; the values for ARC 0-4 are from Yee & 
Lopez-Cruz 1999). The average clustering around the local ULIRGs therefore corresponds 
to an environment similar to the field. A large scatter is seen in our data: although most 
objects are consistent with no galaxy enhancement, a few objects apparently lie in clusters 
of Abell classes and 1. 

Before discussing the results any further, it is important to verify that our analysis of 
the cropped (FOV < 1 Mpc) images does not introduce any bias when compared with the 
results from the uncropped (FOV > 1 Mpc) images. The average (median) B gc -value for 
the 44 objects with uncropped images is 4 ± 121 Mpc L77 (-5 Mpc 1 - 77 ); i.e. slightly smaller 
than the values found for the entire sample. Statistical tests give mixed results regarding the 
significance of this discrepancy (Table 2). The results from a two-sided K-S (Kolmogorow- 
Smirnov) test, a Wilcoxon matched-pairs signed-ranks test, and a Student's t-test on the 
means of the distributions suggest that the distribution of -B 9C -values for the uncropped 
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images is not significantly different from the distribution of -E>c, c -values as a whole, while the 
results from a F-test on the standard deviations of the distributions suggest a significant 
difference. 

We have examined the distributions of -B 9C -values for cropped and uncropped images as 
a function of Galactic latitudes and longitudes. Assuming that the _B 9C -values are unaffected 
by stellar contaminants from our Galaxy, there should be no trend with Galactic latitude or 
longitude. Indeed, the distributions of _B 9C -values with latitude and longitude are consistent 
with being random. Our data therefore confirm the results of Yee & Lopez-Cruz (1999), 
who found that changing the counting radius by a factor of two, both increasing to 1 Mpc 
and decreasing to 0.25 Mpc, did not alter the _B ffC -values significantly. However, given the 
mixed results from the statistical tests, we track the cropped and uncropped data using 
different symbols in the various figures of this paper. We have verified that none of the 
results discussed below are affected if the counting radius is chosen to be 0.25 Mpc or 1 Mpc 
rather than 0.50 Mpc, although quantization errors becomes noticeable when the counting 
radius is 0.25 Mpc due to poorer number statistics. A counting radius of 0.5 Mpc is adopted 
in the rest of this paper to match that of previous studies. 

Next, we explore the possibility of a dependence of the environment richness on ULIRG 
properties. The first parameter we examine is the infrared luminosity (Fig. 1). Statistical 
tests indicate that no trend is present between B gc and Ljf>. This is true for both the 
uncropped data and the entire sample. The same result is found when we examine the 
environment richness as a function of redshift (Fig. 2). Here, however, the redshift range 
covered by our ULIRGs is very narrow (z — 0.1 - 0.22, if we exclude four objects in the 
sample), so this statement is not statistically very significant. Comparisons with the results 
of Blain et al. (2004) and Farrah et al. (2004, 2006) suggest that the environment of high-z 
ULIRGs is richer than that of local ULIRGs. We return to this point in §5. 

In Figure 2, we also distinguish between optical spectral types. We separate our sam- 
ple into Seyfert Is, Seyfert 2s, LINERs, and HII region-like galaxies based on the optical 
classification of Veilleux et al. (1999a). No obvious trends are observed with spectral type, 
but the subdivision of our sample into four subsets necessarily leads to poorer statistics. In 
Figure 3, we plot B gc as a function of log^s/iki); another clear indicator of AGN activity 
[objects with log(/ 2 5//6o) > ~0.7 have "warm", AGN-like IRAS colors]. The lack of trends 
in this figure and Figure 2 indicates that the nature of the dominant energy source in local 
ULIRGs (starburst or AGN) is not influenced by the environment. This result is consistent 
with the ULIRG - QSO evolutionary scenario of Sanders et al. (1988), where the nature of 
the dominant energy source varies with merger phase (starburst in early phases and QSO in 
late phases) but is independent of the environment (as long as the dispersion in velocity of 
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the galaxies within the cluster is not too large to prevent mergers altogether). 

5. Comparison with AGN Samples 

In this section we compare our results with those from published environmental studies 
of AGNs and QSOs. Table 2 summarizes the statistical results of these comparisons and 
Figures 4-9 display the B gc - values from the various samples. Unless otherwise noted in the 
text below, all data sets use the same cosmology. 

5.1. Local Seyferts 

First, we compare our results with those derived on nearby AGNs. De Robertis et al. 
(1998) studied the environments of nearby (z < 0.05) Seyfert galaxies using the exact same 
procedure as the one we use here, so we can directly compare their results with ours. For the 
27 galaxies with z > 0.0045, de Robertis et al. (1998) find (B gc ) = 40 ± 63 Mpc 177 (median 
of 27 Mpc 177 ), consistent with the environment of field galaxies. Recent studies based on the 
SDSS database confirm this result (e.g., Miller et al. 2003; Wake et al. 2004). The average 
environment of local ULIRGs is therefore not dissimilar to that of local Seyferts. However, 
as indicated in Table 2, virtually all statistical tests except perhaps the t-test on the means 
indicate that the two B gc distributions are not drawn from the same parent population. 
Figures 4 and 5 show why that is the case: The distribution of £? 9C -values among ULIRGs 
is distinctly broader than that of the Seyferts. This slight difference is also seen in Figure 6, 
where we display the distribution of local Seyferts and local ULIRGs as a function of Abell 
richness classes. 

As discussed in §4, Seyfert-like ULIRGs do not reside in distinctly poorer or richer envi- 
ronments than non-Seyfert ULIRGs, so the broader scatter in ULIRG environments cannot 
be attributed to the broad range of AGN activity level within the ULIRG population. We 
note that typical error bars for the de Robertis et al. (1998) sample is ~ 100 Mpc L7T , while 
it is ~ 150 Mpc 177 for the ULIRG sample. (The difference is due to the redshift difference 
between the samples - the ULIRG data requires counting to a fainter magnitude, which 
introduces larger uncertainties from background counts.) But the full ULIRG sample distri- 
bution is about 3 times broader than the Seyfert distribution - so, the broader distribution 
of the ULIRG sample cannot be fully explained by the larger error bars. 
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5.2. Local QSOs 

Next, we compare our results with those derived on the nearby (z ~ 0.2) QSOs by Yee 
& Green (1984) and McLure & Dunlop (2001). The measurements of Yee & Green (1984) 
can be directly compared with our ULIRG results since their results were derived using the 
same method and parameters as that of the present study. McLure & Dunlop also apply the 
same formalism to calculate B gc . However, they use a different analysis package to identify 
and classify the objects in the field and carry out the photometry. Their use of HST WFPC2 
data also limits their survey area to only ~ 200 kpc around the QSOs, smaller than even our 
cropped data. These possible caveats should be kept in mind when comparing their results 
with ours. 

Yee & Green (1984) get (B gc ) = 157 ± 208 and a median of 134 Mpc 177 for 34 QSOs 
from the Palomar-Green sample (Schmidt & Green 1983), while McLure & Dunlop (2001) 
derive an average (median) B gc of 365 ± 404 Mpc L77 (241 Mpc L77 ) for a set of 44 radio-quiet 
and radio-loud QSOs and radio galaxies. If we limit our discussion to the QSOs in McLure 
& Dunlop sample (21 radio-quiet QSOs and 13 radio-loud QSOs), the average (median) B gc 
becomes 304 ± 350 Mpc L77 (218 Mpc L77 ). The average environment of the QSOs in both 
studies is therefore slightly richer than that of local ULIRGs. The B gc distributions of the 
two sets of local QSOs (particularly that of the McLure & Dunlop sample; see Fig. 5) show 
a distinct tail at high £? sc -values which is not apparent in the ULIRG distribution. 

A quantitative analysis generally confirms that the B gc distributions of Yee & Green 
and B gc distributions for the radio-quiet and radio-loud QSOs from McLure & Dunlop are 
statistically different from that of the local ULIRGs (Table 2). However, note that the 
Wilcoxon test suggests that the difference is barely significant. Indeed, Figures 5, 7, and 
8 show that there is considerable overlap between the B gc distributions of 1-Jy ULIRGs 
and \ow-z QSOs, particularly the PG QSOs. This result is consistent with the idea that 
some, but perhaps not all, of these QSOs may have formed through a IR-luminous phase 
like that observed at low redshift in the 1-Jy ULIRGs. A more physically meaningful test of 
this scenario would be to compare the environment of local QSOs with the environment of 
z > 0.5 ULIRGs to take into account the finite duration of the ULIRG - QSO evolutionary 
sequence. The recent environmental studies of distant ULIRGs by Blain et al. (2004) and 
Farrah et al. (2004, 2006) indeed point to slightly richer environments, which more strongly 
resemble the environments of the QSOs from McLure & Dunlop. 

A posteriori, the distinct high-_B 9C tail in the distribution of the QSOs of McLure & 
Dunlop (2001) is not unexpected given the host properties of these particular QSOs: ~ 4-5 
times larger host sizes and luminosities relative to the 1-Jy ULIRGs (Dunlop et al. 2003; 
Veilleux et al. 2002, 2006). More luminous hosts live in richer environments on average 
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than hosts of lower luminosity. As pointed out by Veilleux et al. (2006) and Dasyra et al. 
(2006c), the hosts of the QSOs from the Palomar-Green sample (these QSOS are less radio 
and X-ray luminous than the QSOs of McLure & Dunlop 2001) are a better match in host 
size and luminosity to the local ULIRGs. This may explain the generally better (although 
not perfect) agreement between the environments of PG QSOs and 1-Jy ULIRGs. 

5.3. Intermediate-Redshift QSOs 

For the sake of completeness, we display in Figures 5, 6, and 9 the results from our 
study of local ULIRGs alongside the results presented by Ellingson et al. (1991), Wold et al. 
(2001), and Barr et al. (2003) for 63 radio-quiet and radio-loud QSOs at 0.3 < z < 0.6, 20 
radio-quiet QSOs at 0.5 < z < 0.8, and 20 radio-loud QSOs at 0.6 < z < 1.1, respectively. 
All three groups use the same basic method outlined in §3 to calculate the spatial correlation 
amplitude, and all groups assume the same value for H . However, Ellingson et al. (1991) 
assume q = 0.02 instead of 0.5 (Q m = 0.04 instead of 1, if fi> = 0). There is no simple way to 
scale the -B 9C -values for different cosmological models (other than H ) since its computation 
is rather complicated (§3), so Figures 5, 6, and 9 show the B gc - values corrected for the 
different H but not the different Q m . Once again, we see considerable overlap between the 
various distributions, but the statistical analysis formally rules out that they come from the 
same parent population (Table 2). The amount of overlap in B gc - values is quite remarkable 
given the difference in redshifts between the various samples. These results further support 
a connection between ULIRGs and some QSOs. 

6. CONCLUSIONS 

We have derived the spatial cluster-galaxy correlation amplitude, B gc , for 76 z < 0.3 
ULIRGs from the 1-Jy sample and compared our results with those in the literature on 
z < 0.05 AGNs, z ~ 0.2 QSOs, and 0.3 < z < 1 QSOs. The main results are as follows: 

1. Local ULIRGs live in environments which are similar on average to that of field 
galaxies. However, there are a few exceptions: some objects apparently lie in clusters of 
Abell classes and 1. 

2. The infrared luminosity, optical spectral type, and IRAS 25-to-60 /iin flux ratios of 
ULIRGs show no dependence with environment. 

3. The ULIRG environment does not vary systematically over the redshift range covered 
by our sample (mostly 0.1 < z < 0.22). 
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4. There is a lot of overlap between the B gc distribution of local ULIRGs and those 
of local Seyferts, local QSOs, and intermediate- z QSOs. However, quantitative statistical 
comparisons show that the various B gc distributions are not drawn from the same parent 
population. The average environment of ULIRGs appears to be intermediate between that 
of local Seyferts and local QSOs. Local ULIRGs show a broader range of environments than 
local Seyferts, which are exclusively found in the field. The B gc distribution of QSOs show a 
distinct tail at high values which is not seen among local ULIRGs. This slight environmen- 
tal discrepancy between local QSOs and ULIRGs is not unexpected: recent morphological 
studies have found that some of the more radio and X-ray luminous local QSOs used in this 
comparison have more luminous and massive hosts than local ULIRGs. A better match in 
host and environmental properties is seen when the comparison is made with the PG QSOs. 

5. Overall, the results of this study suggest that ULIRGs can be a phase in the lives 
of all types of AGNs and QSOs, but not all moderate-luminosity QSOs may have gone 
through a ULIRG phase. Published studies of the environments of more distant ULIRGs, 
perhaps the actual predecessors of the local QSOs we see today, provide further support for 
an evolutionary connection between ULIRGs and QSOs. 
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Table 1. Galaxies Properties 
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Table 1 — Continued 
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Col (l): Name from the IRAS Faint Source Database. The prefix Z indicates the two objects not in 
the Faint Source Catalog. 

Col. (2): Galactic longitude. 

Col. (3): Galactic latitude. 

Col. (4): Rcdshift from Kim & Sanders (1998). 

Col. (5): Logarithm of the infrared (8—1000 (im) luminosity in units of solar luminosity computed 
using the flux in all four IRAS bands following the prescription of Kim Sz Sanders (1998). 

Col. (6): Optical spectral type from Veilleux ct al. (1999a). 

Col. (7): //MS 25-to-60 fim flux ratio. 

Col. (8): Environment richness parameter computed using PPP program, as described in §3 of this 
paper, in Mpc 1 ' 77 . 

Col. (9): Onc-sigma uncertainty on B gc in Mpc 1 ' 77 . 

Col. (10): Field size in kpc. 



Table 2. Comparisons with AGN and QSO Environmental Studies 
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Col. 


(5): 


Col. 


(6): 



1.77 
1.77 



Sample set used for statistical comparison. 
Number of objects in the sample. 
Mean redshift of sample. 

Mean £? ffC -value and 1-a scatter around the mean of sample in Mpc 
Root-mean square uncertainty on the mean of the S gc -values in Mp 
Median _B gc -value of sample in Mpc 1 ' 77 . 

Cols. (7) and (8): Results from two-sided Kolmogorow-Smirnov test. Entries in col. (7) refer to comparison with the set of 44 ULIRGS that have an image size greater than 1 Mpc X 
1 Mpc, while the entries in col. (8) refer to comparison with the entire set of 76 ULIRG images, regardless of field size. 

Cols. (9) and (10): Results from Wilcoxon matched-pairs signcd-ranks test. Entries in col. (9) refer to comparison with the set of 44 ULIRGS that have an image size greater than 1 
Mpc X 1 Mpc, while the entries in col. (10) refer to comparison with the entire set of 76 ULIRG images, regardless of field size. 

Cols. (11) and (12): Results from Student's t-test on the means of the distributions. Entries in col. (11) refer to comparison with the set of 44 ULIRGS that have an image size 
greater than 1 Mpc X 1 Mpc, while the entries in col. (12) refer to comparison with the entire set of 76 ULIRG images, regardless of field size. 

Cols. (13) and (14): Results from F-test on the standard deviations of the distributions. Entries in col. (13) refer to comparison with the set of 44 ULIRGS that have an image size 
greater than 1 Mpc X 1 Mpc, while the entries in col. (14) refer to comparison with the entire set of 76 ULIRG images, regardless of field size. 



a Thcsc entries refer to the 44 ULIRGS that have an image size greater than 1 Mpc X 1 Mpc. 
^These entries refer to the entire set of 76 ULIR.G images, regardless of field size. 
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Fig. 1. — The environment richness parameter versus the infrared luminosity for local 
ULIRGs. Images which were cropped smaller than 1 Mpc x 1 Mpc are noted by open 
circles, while the non-cropped images are shown as filled circles. The horizontal dashed line 
at B gc = 67.5 Mpc 1 - 77 indicates the average value for typical field galaxies. The range of 
environment richness parameters for the Abell richness classes are marked, following the 
definitions of Yee & Lopez-Cruz (1999). No systematic trend is visible between environment 
richness and infrared luminosity. 

Fig. 2. — The environment richness parameter versus the redshift for local ULIRGs. The 
symbols reflect the optical spectral types of the ULIRGs, as listed in Veilleux et al. (1999a): 
circles are Seyfert 1 galaxies, triangles are Seyfert 2 galaxies, squares are LINERs, and stars 
are HII region-like galaxies. Open and filled symbols stand for cropped and uncropped 
images, respectively. The meaning of the horizontal lines is the same as that in Fig. 1. 
There are no statistically significant trends between environment richness and redshift or 
optical spectral type. 

Fig. 3. — The environment richness parameter versus the logarithm of the IRAS 25-to-60 
jum flux ratio, log(/2s//6o), f° r local ULIRGs. ULIRGs with log(f 25 /f 60 ) > —0.7 are "warm" 
AGN-like systems. The meaning of the horizontal lines and symbols is the same as that in 
Fig. 1. No systematic trend is visible between environment richness and the 25-to-60 /im 
flux ratio. 
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Fig. 4. — Comparison of the environment richness parameters for the local ULIRGs with 
those of z < 0.05 Seyfert galaxies from de Robertis et al. (1998). The meaning of the 
horizontal lines and filled and open circles is the same as that in Fig. 1. Pentagons are the 
data from de Robertis et al. The B gc distribution of local ULIRGs is distinctly broader than 
that of nearby Seyfert s. 

Fig. 5. — Histograms showing the distributions of environment richness parameters for: (a) 
local ULIRGs from this paper (entire sample); (b) local ULIRGs from this paper (uncropped 
data only); (c) z < 0.05 Seyfert galaxies from de Robertis et al. (1998); (d) z 0.2 
PG QSOs from Yee & Green (1984); (e) z « 0.2 QSOs and radio galaxies from Dunlop & 
McLure (2001); (/) z « 0.2 QSOs from Dunlop k McLure (2001); (g) z « 0.2 radio-quiet 
QSOs from Dunlop & McLure (2001); (h) 0.3 < z < 0.6 radio-loud and radio-quiet QSOs 
from Ellingson et al. (1991); {%) 0.5 < z < 0.8 radio-quiet QSOs from Wold et al. (2001; 
model #2 of the background galaxies); and (j) 0.6 < z < 1.1 radio-loud QSOs from Barr et 
al. (2003). The results of statistical comparisons between these various data sets are listed 
in Table 2. None of these data sets appears to be drawn from the same parent population as 
the local ULIRGs, although considerable overlap in the values of the environmental richness 
parameters is seen between the various samples, particularly the local ULIRGs (this paper), 
local Seyferts (de Robertis et al. 1998) and PG QSOs (Yee & Green 1984). 

Fig. 6. — Pie-chart diagrams showing the distributions of environment richness parameters 
typical of field galaxies and clusters of Abell richness classes 0-4 for the eleven different 
samples considered in this paper. See caption to Fig. 5 for a description of these samples. 
The results of statistical comparisons between these various data sets are listed in Table 
2. None of these data sets appears to be drawn from the same parent population as the 
local ULIRGs, although considerable overlap in the values of the environmental richness 
parameters is seen between the various samples, particularly the local ULIRGs (this paper), 
local Seyferts (de Robertis et al. 1998) and PG QSOs (Yee & Green 1984). 

Fig. 7. — Comparison of the environment richness parameter for the local ULIRGs with 
the z 0.2 PG QSOs of Yee & Green (1984). The meaning of the horizontal lines and 
open and filled circles is the same as that in Fig. 1. There is considerable overlap in the 
B gc distributions of local ULIRGs and PG QSOs, although a statistical analysis between 
these two sets of objects generally indicates that they are not drawn from the same parent 
population. 
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Fig. 8. — Comparison of the environment richness parameter for the local ULIRGs with the 
z w 0.2 radio-quiet and radio-loud QSOs of McLure & Dunlop (2001; the radio galaxies are 
not shown). The meaning of the horizontal lines and open and filled circles is the same as 
that in Fig. 1. The environment of these QSOs is distinctly richer on average to that of the 
local ULIRGs, as confirmed in general by a more rigorous statistical analysis. 



Fig. 9. — Comparison of the environment richness parameter for the local ULIRGs with the 
0.3 < z < 0.6 radio-loud and radio-quiet QSOs of Ellingson et al. (1991), the 0.5 < z < 0.8 
radio-quiet QSOs of Wold et al. (2001), and the 0.6 < z < 1.1 radio-loud QSOs of Barr et 
al. (2003). The meaning of the horizontal lines and open and filled circles is the same as 
that in Fig. 1. The B gc distributions of these QSOS overlap considerably with that of the 
local ULIRGs, despite the significant difference in redshifts. 
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